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Abstract 

Percolation and jamming phenomena are investigated for random sequen- 
tial deposition of rectangular needles on d = 2 square lattices. Associated 
thresholds p^^^'^ and p^""^ are determined for various needle sizes. Their ratios 
Pc^^^/Pc'^'" are found to be a constant 0.62 ±0.01 for all sizes. In addition the 
ratio of jamming thresholds for respectively square blocks and needles is also 
found to be a constant 0.79 ± 0.01. These constants exhibit some universal 
connexion in the geometry of jamming and percolation for both anisotropic 
shapes (needles versus square lattices) and isotropic shapes (square blocks on 
square lattices). A universal empirical law is proposed for all three thresholds 
as a function of a. 
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Percolation phenomena are generic in the study of disordered media hke 
porous materials, alloys, ecosystems, etc... A percolation transition is 
based on calculating the probability of occurrence of an infinite connectivity 
between elements of the random medium as a function of the fraction p of 
constitutive elements. The usually concerned object is the percolating cluster 
connecting distant borders of the medium. At the critical point p^'^^^, this 
cluster is very tortuous and is a fractal object. Even though some advances 
have been made in the understanding of percolation, numerous questions 
remains unanswered today 0. Among others, the exact values of many 
thresholds are still not known though some powerful conjectures have been 
proposed for a rather large class of lattices . 

Beside connectivity properties, another fundamental question related to 
disordered structures is the construction of randomly packed structures. 
Among other models, the random deposition also called Random Sequential 
Adsorption (RSA), leads to a critical phenomenon. Above a critical fraction 
pjam q£ occupied sites, the medium becomes fully saturated. No more object 
can be inserted in it [Q. At this stage the packing density of (identical) 
objects filling nearly the whole medium is at maximum. 

In this paper, we investigate both percolation and jamming phenomena 
of non-overlapping anisotropic objects (needles) on strictly d = 2 square lat- 
tices, needles having commensurate scales with respect to the lattice spacing. 
The study of packed anisotropic objects is relevant for geometrical-physical 
properties of granular media like the electrical properties of metallic nee- 
dles P or the compaction properties of granules 0]. We will see that both 
phenomena (percolation and jamming) become highly related to each other 
independently of the aspect ratio a of the needles. 

Earlier studies of needles within the percolation framework have been 
investigated in the overlapping case only and have been introduced for e.g. 
mimicking the formation of microfractures in a brittle material 0. The per- 
colation of non-overlapping needles has never been studied to our knowledge. 
The RSA of needles has only been considered for specific aspect ratio values, 
e.g. a = 2 on d = 2 lattices and for off- lattice cases It seems that 
the RSA of needles on discrete lattices has not been considered up to now. 
The RSA of (isotropic) blocks on a square lattice has been numerically in- 
vestigated by Nakamura [|1T| and thereafter confirmed in |T2|. We will also 
compare our work to Nakamura' s result in order to emphasize the effect of 
anisotropy. Also, an unexpected needle-block relationship will be emphasized 
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in the following! 

On the square lattice, touching needles are needles which have at least 
one cell side in common. Figure 1 illustrates both percolation and jamming 
clusters for the particular aspect ratio a = 4. The figure presents a 16 x 16 
square lattice at •pP^^^'^ and p^"^™- for the case of 4 x 1 needles, and a jammed 
phase for 4x4 blocks. In the case of needles, one observes that both tran- 
sitions take place at quite different values of the fraction p of occupied sites. 
As we will see below, the percolation threshold is "low" implying that the 
unoccupied sites (holes) can form open and large pores. At the percolation 
threshold, the needle structure is highly heterogeneous. When the fraction 
of needles increases, the large pores are filled up to the jamming situation. 
Near the jamming transition, some short range order appears. Indeed, a 
close packing of needles oriented in the same direction is observed. Close 
(oriented) packed needles seem to form blocks of size a x a! In order to em- 
phasize this ordering, a typical jamming phase is presented in Figure 2 for 
L = 400 and for a = 20. Horizontal and vertical orientations of needles are 
represented in different grey levels. The different colors put into evidence the 
short range ordering. More precisely, some domains of horizontal or vertical 
needles are seen. 

In the third case illustrated in Figure 1, 4 x 4 blocks are seen in the 
jamming phase. Large spanning clusters of connected blocks are not formed. 
One should note that percolation does not occur below the jamming threshold 
in the case of blocks. However, percolation always takes place before jamming 
in the case of needles. This makes a big difference between the deposition of 
anisotropic and isotropic non-overlapping objects on square lattices. 

In order to determine the thresholds for both phenomena, numerical sim- 
ulations have been performed on d = 2 square lattices containing up to 
2000 X 2000 sites. Starting from an empty lattice, needles are added sequen- 
tially such that p increases linearly. The percolation and jamming phenom- 
ena are checked until they are found. The first percolation threshold, i.e. 
when two distant borders are connected, is considered here, in contrast to 
the second percolation threshold which is obtained when the four borders are 
connected |jl3[- The probability Pperc to find a percolating cluster and the 



probability Pjam to find a jamming phase are fitted by the error function 
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where pc is the critical point and A is the width of the transition. Our 
assumption that the distribution of critical points is a Gaussian is sufficient 
from a practical point of view, though not claimed to be exact 0. One should 
also note that the jamming transition is always sharper than the percolation 
transition, i.e. /S.jam < ^perc- The length A vanishes as a power of the 
system size L. One has 

A ~ L-^/" (2) 

where ~ means asymptotic proportionality. This allows for the measure of 
the exponent u for the correlation length ^ which diverges at the critical point 
as ^ ~ \Pc — p\~'^ ■ The above relationship (2) means that the transition is 
sharper for larger lattice sizes. As expected, the exponent value has been 
found to be i^perc = 1-35 ib 0.02 compatible with the I'perc = 4/3 value known 
for d = 2 percolation with isotropic particles. Another exponent value is 
found for the jamming transition: Ujam = 1-0 ±0.1, a value also reported 



in the work of Nakamura |Tl| about the jamming of blocks. Both values 
Uperc and Ujam havc been found to be independent of the aspect ratio a of 
the needles. It seems at first that both critical phenomena are independent 
of each other since critical exponents and thresholds are different. However, 
we will see herebelow that a relationship between percolation and jamming 
thresholds nevertheless exists! 

The threshold values are slightly dependent on the system size L. The 
following dependence 

p,(oo)-p,(L)~L-i/^ (3) 

has been established for percolation IQ. The latter relationship allows us to 
extrapolate the threshold for an infinite system L —>■ oo. Table I summa- 
rizes some of the extrapolated values of pc for both percolation and jamming 
phenomena for typical values of a. When the aspect ratio a of the needles 
increases, both thresholds expectedly decrease. Figure 3 presents both p^'^^'^ 
and pi"'"^ as a function of a. Curves are guides for the eye. 

In order to search for a possible relationship between both thresholds if 
any, we have calculated the ratio p^'^^'^ / pi"-"^ (see the third column of Table 
I). Surprisingly, this latter ratio is found to be a constant 0.62 ± 0.01 what- 
ever the aspect ratio value a! Small deviations from the 0.62 value are only 
observed for large a values; those deviations are certainly due to the parti- 
cle/lattice finite size. The universal ratio is unexpected and suggests that 
both critical phenomena are intimately related, although we have found that 
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the critical exponent values are different. In fact, this result indicates that 
the percolation cluster should be a fundamental cluster (like a skeleton) for 
the jamming phase! 

It should be noted that the scaling properties of the jamming phase are 
not similar to those of dense systems otherwise the threshold p^"*" would 
be independent of a. In the case of the jamming transition, a power law 
behavior 

P^r - a-' (4) 

with an exponent 6 ~ 0.20 has been proposed for off-lattice RSA of rectangles 
1^ . These authors [[10[ have suggested that the power law (4) indicates that 



the network of needles is fractal with a dimension 2 — 5 at p^"™ at scales below 
a. Moreover, systematic deviations from Eq.(4) appear for large a values. It 
seems that both thresholds are decaying in a slower fashion than a power of a 
(not shown in Figure 3). An empirical law will be proposed in the following. 

Consider now the Nakamura work. The ratio between the thresholds 
pjam blocks and needles seems also to be a universal constant (see the 
last column of Table I)! This remarkable result suggests that there is also a 
relationship between needles and blocks for the jamming phase. Of course, 
we have seen above that close packed needles form clusters which can be seen 
as blocks of size a x a. 

Let us interpret our results by considering a coarse grained view of the 
disordered systems. In a coarse grained view, i.e. at a scale larger than a, 
the notion of anisotropy should disappear. Needles of size a are replaced by 
blocks of size a x a. The threshold is then expected to be that of classical 
square objects. Since the notion of anisotropy disappears at larger scales, 
the fraction of occupied sites p' in a coarse grained view is larger than the 
true fraction of occupied sites p. We propose the equation 



p = p' 



7 



a — 1 



(5) 



where 7 is a constant. The factor in the r.h.s. represents the fraction of 
free sites remaining if two perpendicular needles are touching in any a x a 
"supersite" . This term represents the loss of information when one looks for 
the connectivity of the packing in a coarse grained view. This term represents 
also the jamming phase since no more needles can be added in a "supersite" 
which contains already 2 perpendicular touching needles. The parameter 7 
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should depend only on geometrical aspects and should be independent of a. 
Using the relationship (5), one expects that both thresholds scale as 



where C is a constant. The law (6) fits the data and is illustrated in Figure 4. 
The agreement is quite remarkable. The coefficient 7 is found: 7 = 0.31±0.01 
in all cases, providing a constant ratio of thresholds. One should note that 
the Eq.(6) is not valid for a = 1 which is a particular (isotropic needles!) 
point. 

In the Nakamura case (blocks of size a) , one has a loss of information due 
to block edges. The empirical law of Eq. (6) is found to hold also in that case! 
This good agreement suggests that our physical arguments are appropriate 
to the study of packed anisotropic objects. 

In summary, we have investigated two phenomena, i.e. percolation and 
jamming, which have been independently studied up to now. We have found 
that they are closely related in the case of anisotropic objects (needles). We 
have interpreted this effect assuming a loss of information in a coarse grained 
view. Fundamental laws have been obtained theoretically. Finite size and 
shape effects seem to have been captured in a scaling law of (1 — 1/a)^ for 
needles. 
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a 




needles 

pjam 


perc 
Pc 

Pc 


blocks 

pjam 


d""" (blocks) 

pl"'"^ (needles) 


2 


0.562 


0.906 


0.621 


0.749 


0.827 


3 


0.528 


0.847 


0.623 


0.681 


0.804 


4 


0.504 


0.811 


0.622 


0.646 


0.797 


5 


0.490 


0.787 


0.623 


0.628 


0.798 


6 


0.479 


0.770 


0.622 


0.620 


0.805 


8 


0.474 


0.757 


0.626 


0.603 


0.797 


10 


0.467 


0.741 


0.630 


0.593 


0.800 



Table I — List of some thresholds Pc for both percolation and jam- 
ming phenomena for typical values of the needle aspect ratio a. The ratio 
pperc ipjam jg ^^so given in the fourth column. The Nakamura thresholds for 
blocks in the jamming phase are given in the fifth column. In that case, a 
represents the block width. The last column gives the ratios for blocks and 
needles thresholds p^""*. 
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Figure captions 



Figure 1 — (a) Percolation cluster of non-overlapping needles of aspect 
ration a — 4. (b) The largest cluster at the jamming transition. 

Figure 2 — Jamming phase for a — 20 needles. The lattice size is 500 x 500. 
Vertical and horizontal orientations are represented by different grey levels. 

Figure 3 — The thresholds p^^^"^ and p^"™ as a function of the needle aspect 
ratio a. The nakamura results are also illustrated. The curves are only guides 
for the eye. 

Figure 4 — The thresholds p^^^^ and p^""* as a function of (a — l)^/a^. The 
Nakamura results are also illustrated. The straight lines are fits using Eq.(6). 



8 



References 

[1] D.Stauffer and A.Aharony, Introduction to Percolation Theory, (Taylor 
& Francis, London, 1994) 2nd printing 

[2] D.Stauffer, Physica A 242, 1 (1997) 

[3] S.Galam and A.Mauger, Phys. Rev. E 53, 2177 (1996); S.C. van der 
Marck, Int. J. Mod. Phys. C 9, 529 (1998) 

[4] J.W.Evans, Rev. Mod. Phys. 65, 1281 (1993) 

[5] E.Pirard, private communication 

[6] E.Caglioti, V.Loreto, H.J.Herrmann and M.Nicodemi, Phys. Rev. Lett. 
79, 1575 (1997) 

[7] P.C.Robinson, J. Phys. A 17, 2823 (1984) 

[8] R.S.Nord and J.W.Evans, J. Chem. Phys. 82, 2795 (1985) 

[9] R.D.Vigil and R.M.Ziff, J. Chem. Phys. 91, 2599 (1989) 

[10] R.M.Ziff and R.D.Vigil, J. Phys. A 23, 5103 (1990) 

[11] M.Nakamura, J. Phys. A 19, 2345 (1986) 

[12] V.Privman, J.S.Wang and P.Nielaba, Phys. Rev. B 43, 3366 (1991) 
[13] M.Ausloos and J.M.Kowalski, Phys. Rev. B 45, 12380 (1992) 



9 



